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Abstract. The presence of a primordial magnetic field would have induced resonant 
conversions between photons and axion-like particles (ALPs) during the thermal 
history of the Universe. These conversions would have distorted the blackbody 
spectrum of the cosmic microwave background (CMB). In this context, we derive 
bounds on the photon- ALP resonant conversions using the high precision CMB spectral 
data collected by the FIRAS instrument on board of the Cosmic Background Explorer. 
We obtain upper limits on the product of the photon- ALP coupling constant g times 
the magnetic field strength B down to gB <10~^^ GeV~^nG for ALP masses below 
the eV scale. 
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1. Introduction 

The possible existence of a primordial magnetic field of cosmological origin has been the 
subject of an intense investigation during the last few decades. Despite these efforts, 
no astrophysical evidence has been reported so far concerning magnetic fields over 
cosmological scales, and only upper limits are reported. When scaling the original bound 
from the Faraday effect of distant radio sources [1,2] to the now much better known 
baryon density measured by the Wilkinson Microwave Anisotropy Probe (WMAP) [3], 
one has B <2.8 x 10"''(/c/Mpc)~^/^ G, coherent on a scale Ic ^ 1 Mpc [4]. A recent 
analysis of the WMAP 5-year data on the Faraday rotation of the linear polarization of 
the cosmic microwave background (CMB) gives comparable upper limits ranging from 
6 X to 2 X 10-6 G [5]. 

The presence of a primordial magnetic field would inevitably produce resonant 
conversions between photons and axion-like particles (ALPs) in the Early Universe. 
This effect was pointed out at first in a work by Yanagida and Yoshimura [6]. They 
realized that such a conversion could produce a sizable distortion in the CMB spectrum. 
Since at that time data indicated a distortion in the Wien region of the CMB, these 
authors speculated that a resonant axion-photon conversion could be an intriguing 
explanation for such an effect. Nowadays, the blackbody nature of the CMB spectrum 
has been measured with a precision better than 1 part in 10^ by the Far Infrared 
Absolute Spectrophotometer (FIRAS) on board of the Cosmic Background Explorer 
(COBE) [7, 8]. The astonishing accuracy of such a measurement now allows to constrain 
exotic scenarios which would deplete the CMB spectrum, such as hidden photons [9, 10], 
axions [11], radiative neutrino decays [12] or millicharged particles [13]. Recent bounds 
on photon- ALP mixing using the precision CMB data [11] have focused only on non- 
resonant conversions in the late Universe. However, for values of ALP masses allowing 
resonant conversions during the expansion of the Universe one expects stronger bounds 
due to this resonant enhancement. For this reason, in the present paper we determine 
the bounds from resonant photon-ALP conversions in light of the COBE precision 
CMB data. Moreover, we make use of an accurate description of the evolution of the 
cosmological plasma in which the resonant conversions arise, following the line of our 
recent study on resonant photon-hidden photon conversions in the Early Universe [10]. 

Often, strong bounds on the photon-ALP coupling arise from stellar evolution [14]. 
Our bounds are difficult to compare with these other constraints since they strongly 
depend on the strength of the primordial fields. If fields near the experimental bounds 
~ 100 nC are realized, our bounds can be much stronger than stellar evolution 
constraints, but they could also be weaker if the primordial fields were much smaller. 
However, in general, our bounds are complementary to stellar evolution ones, since some 
ALP models [15-19] predict a suppression of the photon-ALP coupling in dense stellar 
interiors, which would not happen in the relatively diluted primordial plasma. For this 
reason, we think worthwhile to explore which independent constraints are achievable 
from cosmological arguments. 
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The plan of our work is as follows. In Section 2 we review the mechanism of mixing 
between photons and axion-like particles and we present our analytical prescription to 
calculate the resonant photon-ALP conversion probability in the expanding Universe. 
In Section 3 we describe our simplified model for the effective photon mass induced 
by the primordial plasma. In Section 4 we discuss the effect of the random primordial 
magnetic fields on the resonant conversions and we characterize the averaging of the 
conversion probability over the sky and over the photon polarization, relevant to study 
ALP effects on the CMB monopole spectrum. In Section 5 we describe the constraints 
coming from spectral CMB distortions for ALP masses undergoing resonant conversions 
after the recombination epoch. In Section 6 we extend our limits to pre-recombination 
resonances, using the experimental limits on the chemical potential fi of the CMB 
spectrum and the agreement between Big Bang Nucleosynthesis (BBN) and the CMB 
on the effective number of additional relativistic species at decoupling. Finally, in 
Section 7 we comment on the complementarity of our bound with other astrophysical 
and experimental constraints and present our conclusions. 



2. Photon mixing with axion-like particles 

2.1. Axion-like Particles 

Axion-like particles (ALPs) are scalar or pseudoscalar bosons (f) that couple to two 
photons with one of the following interaction Lagrangians [20] 

^pseudoscalar = — ^Q-F^ivF^" (j) = Q-B ■ E0 , (1) 

Csc.i.r = \g+F,,F^''4>=^g+{B'-E^)4> , (2) 

where F'^'^ = e'^'^"^F„^/2 is the dual of the field strength tensor and E and B are 
the electric and magnetic fields, respectively. In the presence of a constant external 
magnetic field one can decompose the electromagnetic field into an external component 
and the dynamical part representing photons as B — B''^* + V x A. The terms in the 
Lagrangian containing B^^* act as mass mixing terms between the axion-like particle 
and the photon field A^^ = {Aq; A). Choosing the radiation gauge V • A = the mixing 
terms become evident, 

^pseudoscalar = fi'-B*'''* ■ A + ... , (3) 
^scalar = ^?+B°^* ■ ( V X A)0 + . . . . (4) 

In our study we will consider the mixing of photons with these particles. If we write A 
as a plane wave with frequency u and wave vector k we can explicitly evaluate its time 
and spatial derivatives. Reabsorbing a factor of i into A and using uj ~ |k| we find that 
the mixing part in the Lagrangian reads 

£-iar =g+BTu;A^<p, (6) 
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where is the component of the external magnetic field perpendicular to the 
propagation direction of photons (k) and yl||,y4x are respectively the components of 
A parallel and perpendicular to that component. Note that the component of B*^^* 
parallel to k does not contribute to the mixing and also that pseudoscalar fields mix 
with the parallel component A\\, while scalars mix with A±. In the following, we denote 
with g the ALP-photon coupling constant and consider only the photon component, 7, 
that mixes with 0. 

Due to the effective mass-mixing in the external magnetic field, the propagating 
eigenstates in vacuum are now rotated with respect to (7, 0) by an angle 6 given by [21] 

sm^t^ = , , (7) 

/m^ + {2gBujy 

9 

771/ 

cos 2^ = = , (8) 

/m^ + {2gBuy 



where, for simplicity, we have indicated the transverse component of the magnetic field 
with B, and is the ALP mass. This misalignment is well known to produce 7^0 
oscillations with a wavenumber given by [21] 



^ ^ /m^ + i2gBu;)^ ^ Am^ 
~ 2u ~ 2uj ' ^ 

2.2. Medium effects and resonant MSW conversions 

Photon oscillations into ALPs are modified by the refractive properties of the medium. 
In the primordial plasma, photons acquire a non-trivial dispersion relation which can be 
parametrized by adding an effective photon mass to the Lagrangian. This is generally 
complex, reflecting the absorption properties of the plasma. However, for most of the 
parameter space studied here it turns out that absorption is negligible. Moreover, even 
when this is not the case our results will show no dependence on it. Therefore, we have 
chosen to neglect it in the following exposition and include a few remarks when relevant. 
In this case, the effective mixing angle 9 is related to the vacuum one by [14] 

sin2^~ = ^ , (10) 

sin^2e + {cos2e-if\ 

~ cos2^-f 
cos2^ = ^ ' (11) 

sin2 2^ + (cos2^-0^^ 



where the parameter ^ which measures the signiflcance of the medium effects reads 

=co.2e{'<] . (12) 
lm\ + {2gBuf \m^J 

As the Universe expands, eventually the condition 

= , (13) 
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is satisfied and 6 it /A. When this condition is fulfilled, resonant photon- ALP 
conversions are possible, analogous to the well-known Mikheev-Smirnov-Wolfenstein 
(MSW) effect for neutrino flavor transitions [22-24]. If the photon production and 
detection points are separated by many oscillation lengths from a resonance, the 
oscillation patterns wash out. Thus, the transition probability is given by [25] 

P-y^^ (p-l-) COS 29o COS 29 , (14) 



2 V 2, 

where we have assumed the mixing angle 6o at the detection in vacuum, 6 is the effective 
mixing angle at the production point (considered to be at high density) and p is the 
level crossing probability. This latter takes into account the deviation from adiabaticity 
of photon- ALP conversions in the resonance region. In particular, one has p = for 
a completely adiabatic transition and p = 1 for an extremely nonadiabatic one. The 
crossing probability p for photon- ALP resonant conversions can be obtained using the 
Landau-Zener expression [10] 

p ~ exp (^—2'iTrk sin^ 6'^^ , (15) 

where k is again the 7—^0 vacuum oscillation wavenumber in Eq. (9), 6r is the vacuum 
mixing angle at the resonance and 

d\nm'i{t) 



dt 



(16) 

t — ^rcs 

is a scale parameter to be evaluated at the location where a resonance occurs. $ 

Note that the mixing angle 6 changes as the Universe evolves because the frequency 
u and the magnetic field B are functions of the redshift. Denoting with uq and Bq 
their values at 2; = 0, the frequency grows as uj = uJo{l ~\- z) while the magnetic field 
evolution is mo del- dependent. In this paper we have decided to focus in the most 
widely studied case of considering magnetic fields frozen into the medium, for which 
B = Bq{1 + z)"^ [2]. It is clear from Eq. (7) that at early times the vacuum mixing 
angle becomes maximal, sin 26' 1. However, the effective mixing angle 6 remains 
small because of the suppression by matter effects. In fact, at early times, before 
recombination, is given by the plasma frequency u"^ which is proportional to the free 
electron density which scales as (1 + z)^. Therefore, at high redshifts {z 00), ^ tends 
indeed to a constant given by 



m. 



2 



e-^ = ^^^^^-2.7xlO^(,,o5o,nG)-M^ , (17) 

where we have used my{z = 0) = 1.59 x 10~^^ eV [10], normalized the photon frequency 
today ujq with respect to the actual CMB temperature Tq = 2.725 K [8] and finally 
defined gio = g/10^^^ GeV~^ and i?o,nG = Bq/I nG as typical values close to their 

t Our problem is different to the widely discussed of neutrino oscillations or the photon-hidden photon 
case studied in [10] since here even the vacuum mixing evolves in time. However, the corrections to 
Eq. (15) that we are neglecting are unimportant, since the evolution of the vacuum mixing is very 
smooth near the resonance, contrarily to the evolution of 9. 
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current experimental limits. For such large value of ^ we can take cos 2^ ~ —1 at the 
production point in Eq. (14). 

Moreover, the FIRAS sensitivity will allow us to bound P^^^^IO"^ which also 
excludes non-resonant vacuum oscillations with sin^ 29o > 10~^. Considering resonant 
conversions we will constrain much smaller mixings so that we can also take cos 26*0 — 1 
in Eq. (14). This then simplifies to 

P^^^ ~ 1 - p . (18) 

The FIRAS bounds thus require a strongly non-adiabatic resonance. Under this 
condition, we can approximate the conversion probability as 

P^_^ ~ 27rrA;sin2^^ . (19) 

The product rk is basically the ratio of the characteristic expansion time of the universe 
and the vacuum oscillation length, and for the region of parameters considered 

here is large. Since FIRAS constrains P^^^ <^ 10"^, this tells us that the vacuum mixing 
angle Or at the resonance has to be extremely small. Therefore, we will work in the 
small mixing regime, for which Am^ 2^ and 

uJo\ /10-^^eV^^ 



sin 2^, ~ 0.92 X 10"^ (1 + z)' igioBo,nG) [yJ [ ^a^ ) ' (^0) 

where in Eq. (7) we have approximated 6,. ~ gBu/m'^, since at the resonance = 
and ^ 2gBuj for the values of ALP masses we are considering. Under, this 
approximation, the expression for the conversion probability is exceedingly simple, i.e. 

q'^B'^Tiruj , , 

ml 



3. Cosmological profile 

In order to calculate the conversion probability we need the profile of the photon effective 
mass along the cosmological line of sight. In this work we use the same prescription 
of our Ref. [10] in which we studied the resonant conversions between photons and 
hidden photons in the expanding Universe. 

The effective mass can be parametrized as 

= cuJ(Xe) X 

where uj^^Xe) ~ 1.59 x 10~^''(1 + zY/'^X]/'^ eV is the plasma frequency with Xe{z) the 
hydrogen ionization fraction as a function of redshift, which we take from [27]. We note 
that the CMB energies probed by FIRAS are much smaller than the first excitation 
energies of the hydrogen (when the temperatures are low enough to allow bound atoms) 
so that we have neglected these effects in the expression of the effective mass. We 
have also ignored the subleading helium contribution. Based on [10], in Fig. 1 we show 
our reference profile for the cosmological history of the effective photon mass as a 



0,0073 I -J 



U \2 /I -X, 



(22) 
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10^' ^ lO' 10^ 10^ 

z 

Figure 1. Effective photon mass as function of rcdsliift. The thin line uses — I, the 
blue, magenta, green and red lines are for co/T = 1, 3, 4, 10, respectively. These values 
are in the frequency range probed by FIRAS. The two sharp dips in the green and red 
lines bound the region where becomes negative. (See [10] for further details) 

function of redshift z. We recall that the evolution of the effective photon mass which 
depends on the history of the ionization fraction Xe{z) and also on the photon energy, 
is extremely complex. Above a temperature T ~ 0.5 eV (redshift z ~ 1100) hydrogen is 
fully ionized. As the Universe temperature decreases, photons cannot ionize hydrogen 
efficiently and electrons and protons slowly combine. This makes the Universe very 
transparent to radiation, indeed releasing the photon bath which we see today as the 
CMB. This is the so-called recombination epoch. Later on {z <, 6) the Universe becomes 
ionized again due to ultraviolet radiation from the first quasars or population III stars. 
These effects produce the non-monotonic behavior in the shape of the effective photon 
mass as a function of redshift seen in Fig. 1. Note also that during recombination the 
effective mass squared can become negative if the photon energy is large enough. 

Along such a non-monotonic mass profile multiple resonances can occur for a given 
ALP mass (provided < 10~^^ eV). If this is the case, Eq. (21) tells us that the most 
adiabatic crossing, quantitatively the more important, is the one that happens at earlier 
times. This is particularly clear when we express all red-shift dependent quantities as 
their values at z = times their redshift dependence and use r oc {3H)~^, with H the 
Hubble constant, to find Pj-kp oc {1 + z)^/H which increases with redshift since H grows 
at most as (1 + 2;)^ during the radiation dominated epoch. 

We note that for photon energies tu/T^S.S there is a small period during the dark 
ages for which the effective mass squared becomes negative. This implies that there 
is one redshift at which = and at which, in principle, ALPs of arbitrarily small 
mass can be resonantly produced. However, due to the fast drop of m^, this crossing 
is extremely non-adiabatic, as we have explicitly checked. For this reason, we will not 
consider this case hereafter. 
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In summary, in the case of multiple crossings the most relevant one happens before 
reionization. As a consequence, for > 10~^^ eV we need to evaluate resonances 
only for 2 ^ 10. This situation simplifies our calculations because after reionization the 
plasma is complicated by the presence of density inhomogeneities [10]. 



4. Magnetic fields at the resonance and probability averages 

In principle, the primordial magnetic field is not known, although it is generally expected 
to show a somewhat turbulent structure [28]. We would like to clarify how it affects 
the resonant photon-ALP conversions. For this reason, we have to compare the typical 
width of a resonance region with the coherence scale of the magnetic field. The half- 
width of the resonance is, according to Eq. (10), S^{t) = sm26j-, which corresponds to 
a length scale 

Tr = r sin 26r (23) 

^ 1.3 X 10-2(1 + ^)3/2 (g. 




^ioSo,nG ^ 1.3 X 10-2 ( ) ( ^ ) ( ^ 1 . (25) 



^ 0.42 ((7ioSo,nG) j [ ^\^^' 1 Mpc , (24) 

where for the first numerical estimate we have used Eq. (20) and r < {3H)~^ ~ 1.4(1 + 
z)-3/2Qpc^ with H the Hubble constant. Eq. (24) results from the fact that for the 
most relevant resonance before reionization (at z>10) one has = m^>10"^^ [(1 + 
z)/10f/^eV. The resonance half-width is thus smaller than the coherence length ~ 
h,o (1 + of the magnetic field for 

_2 f y/3 /To\ ( lc,o ' 

Wo/ VMpc^ 

The boundary of our exclusion bounds will satisfy this constraint (see Fig. 3) so 
that we can consider the magnetic field constant during the resonance to compute it. 
At couplings larger than the ones at the exclusion boundary, the resonance width will 
eventually become larger than the magnetic field coherence length and photons will 
see different magnetic field domains during the resonance. Nevertheless, the transition 
probability increases by increasing the value of g, so that we can also exclude this 
upper region. Therefore, also in such a situation, we will continue to use our simple 
prescription of taking the magnetic field evaluated exactly at the crossing point, defined 
by our Eq. (13), and constant during the resonance§. 

However, the magnetic field direction and strength during the resonance will 
generally be different along different directions in the sky. Moreover, the photon-ALP 
resonant conversions depend on the relative orientation of the photon polarization and 
the magnetic field direction, so that in different magnetic domains different photon 

§ Possible corrections due to the variation of B inside of the resonant region (see, e.g., the discussion 
in [26] in the neutrino case) are not expected to modify our exclusion plot above the boundary. 
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post-recombination 



pre- recombination 



Figure 2. Picture of resonant photon- ALP oscihations of the cosmic microwave 
radiation. The transitions happen during a very short time in the history of 
the universe (the resonant sheU) determined by the ALP mass which here 
is shown confined by dashed hnes. The primordial magnetic field can be 
assumed to be constant during the resonance (the shell thickness) but both its 
strength and orientation depend on the angular position in the sky. Resonant 
transitions happening after recombination deplete CMB photons depending 
on their polarization. If resonances happen before recombination, photon 
trajectories can form any angle with B during the resonance and then be re- 
scattered towards the observer such that information about the polarization is 
lost. In this case, the CMB light from the last scattering surface is effectively 
averaged over the angle ip. 



polarization states play the role of A\\ and A^, see Fig. 2. For a generic photon 
polarization, the B strength entering the conversion probability in Eq. (21) is|| 



where x is the position vector of the resonance region in a particular direction x, e is the 
photon polarization vector (|e| = l,e x x = 0), ipix) is the angle between the magnetic 
field B(x) and the photon propagation direction x and ip the angle between Bt (the 
component of the magnetic field perpendicular to x) and e, see Fig. 2. 

In this paper we will mainly discuss bounds of photon-ALP mixing from the CMB 
monopole. If the resonance happens after recombination, the full-sky and polarization 
averaged photon-ALP conversion probability Eq. (21) is 



II Eq. (26) holds for parity-odd ALPs; when considering parity-even one finds |B(x) x e| = 
|B(x) sin-!/)(x) sin(^| instead. 




(26) 




(27) 
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where we have defined a conveniently normahzed average 

{B) = 3 / / — B sm -wcos 09 

^ ^ i 47r Jo 27r' ' ^ ^ 



^ ^"^^^iBpsin^^ , (28) 



47r 

when we recall that implicitly B and ip depend on x. The convenience of this average 
is clarified in two simple cases: a) the absolute value |B| is constant in the whole sky 
but has a random direction and b) it has a fixed direction. In both cases we find 

(B^) = \Bf . (29) 

When the resonance happens before recombination, the photon trajectories can 
form any angle ip with B during the resonance and then be re-scattered towards 
the observer such that information about the polarization is lost. Since the photon 
trajectories are not straight, the angle is not correlated with the magnetic field 
direction, photon polarization and the direction in the sky. Because of this, one has to 
perform a local average over ip in each magnetic domain before doing the sky-average, 
i.e. 

2 [dipdVt^ . 2 , 2 1|-D|2 (nr,\ 

|B| / — ^ sm COS <y5 = -|B| . (30) 

J 2% An 3 

Assuming that B = B(x) does not vary very much on the distance scales of the horizon 

at last-scattering surface, we shall then define a sky average 

{B') = / ^|BP , (31) 

so that also in the pre-recombination epoch one can use Eq. (27) for the conversion 
probability. Note that if |B| is isotropic, also in this case we obtain (5^) = |Bp. This 
expression neglects in principle decoherence effects from photon scattering absorption 
during the resonance. However as we will see in Sec. 6 these effects are not relevant if 
the conversion probability is small. 



5. Bounds in the post-recombination era 

The CMB spectrum measured by FIRAS fits extremely well to a black-body spectrum 
at a temperature Tq = 2.725 ± 0.002K [8]. The energy range of the CMB spectrum 
measured by FIRAS [7] is 2.84 x 10"^ eV < ujq < 2.65 x 10~^eV, corresponding to 
1.2 < uiq/Tq < 11.3. In that region, the CMB blackbody becomes unprotected to 
distortions below a cosmic temperature ~ 50 eV, which corresponds to a photon plasma 
mass of ~ 10^^ eV. On the other hand, today the average plasma mass for photons 
is as low as 2 X lO"-"^"^ eV. If ALPs exist with a mass between these two values they will 
be produced resonantly and leave their imprint on the CMB as a frequency-dependent 
distortion [see Eq. (27)]. Bounds for higher masses are considered in the following 
section. 

In order to obtain our bound, we have considered the distortion of the overall 
blackbody spectrum. To this end we use the COBE-FIRAS data for the experimentally 
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measured spectrum, corrected for foregrounds [7]. Note that the new cahbration of 
FIRAS [8] is within the old errors and would not change any of our conclusions. The 
= 43 data points ^f^^ at different frequencies Ui are obtained by summing the 
best-fit blackbody spectrum to the residuals reported in Ref. [7]. The errors a^^^ are 
also available. In the presence of photon-ALP conversion, the original intensity of the 
"theoretical blackbody" monopole at temperature T, 

$°(c.,T) = i4[exp(c./T)-l]-^ , (32) 

would be deformed to 

<l>(a;,T,A) = $°(c.^,T)[l-(P^_^)] , (33) 

where (P^^^) is the sky average of the polarization averaged photon-ALP conversion 
probability, defined in Eq. (27). We can then build the reduced chi-squared function 



1 



N 



$r-$(^„T,A) 



2 



(34) 



We minimize this function with respect to T for each point in the parameter space 
A = (m^, (7(5^)^/^), i.e. T is an empirical parameter determined by the xt minimization 
for each A rather than being fixed at the standard value. 

In Fig. 3 we show our exclusion contour. In particular, the region above the 
continuous curve is the excluded region at 95% C.L., i.e. in this region the chance 
probability to experimentally obtain larger values of xt is lower than 5%. 

We shall stress again that the CMB photons from different angles have traversed 
magnetic fields during the resonance that have in principle different strengths and 
directions. This would induce anisotropies and polarisation. The difference between 
a direction parallel and orthogonal to the magnetic field must not exceed the observed 
CMB temperature anisotropy, AT/T < 10"^. However, the bound achievable from this 
anisotropy pattern is expected to be less stringent than the one put from the distortion of 
the overall blackbody spectrum [1 1] . Possible improvement of this bound would require 
a detailed investigation in terms of the multipole expansion of the CMB temperature 
fluctuations. We deserve this task for a future work. 



6. Bounds in the pre-recombination era 

The resonant conversion of photons into ALPs produces an energy-dependent depletion 
of the CMB which can be constrained by FIRAS data. However, if this resonance 
happens before recombination, the CMB is still coupled to the primordial plasma and 
these distortions can be processed by photon-plasma interactions which will tend to 
thermalize the spectrum. Photon scattering and absorption during the resonance can in 
principle affect our result Eq. (21) producing a damping of the photon-ALP conversions. 
Both effects can be introduced in an imaginary contribution to the effective photon mass 
squared. In the limit of strong damping, the flavor relaxation rate becomes independent 
on the details of the scattering processes. When this situation dominates a resonance. 
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taking the results of [9] (computed in the photon-hidden photon case but vahd for a 
general mixing case) and translating them into our problem in the relevant small mixing 
case, we obtain 



where p is the same crossing probability given by Eq. (15). The effects of decoherence 
become manifest in the fact that for adiabatic resonances, p ~ 0, the system reaches 
thermal equilibrium, and photons and ALPs equilibrate their populations which, of 
course, results in a transition probability of 1/2. However, in the non-adiabatic limit, 
this expression gives exactly our Eq. (21), with r = {3H)~^ in this regime, so our results 
do not depend on whether the resonance is vacuum or damping dominated and we can 
use Eq. (21) in both the cases. For the sake of the simplicity, we neglect possible cases 
in which the scales of the damping and of the resonance are of the same order. As 
a summary, we only have to consider the processing of the CMB spectrum after the 
resonance to compare it with the FIRAS data. 

Note that the processing strongly depends on the temperature at which the 
resonance takes place, since it determines if the different photon interactions are effective 
or not. The response of the primordial plasma to distortions of the Planck distribution 
is driven by two fundamental classes of processes: Compton scattering conserves the 
number of photons, whereas double Compton scattering, Bremsstrahlung and their time 
reversal analogues change the number of photons. 

Compton scattering is the fastest process and is responsible for the change of photon 
direction and polarization described already in Sec. 4. Compton scattering can also be 
effective at redistributing the energies of photons, however per se cannot change the 
photon number. It can provide kinetic equilibrium of the photon distribution provided 
it is efficient, i.e.. 



where ciTh is the Thompson's scattering cross section. This condition corresponds to 
masses ^ 1.4 x 10~^ eV. Therefore, below this mass the plasma is so weakly coupled 
that it cannot process the distortion. In such a condition, the treatment given in the 
previous section holds with the only difference being the averaging of photon directions 
during the resonance, as discussed in Sec. 4. The corresponding bound is shown in 



For masses ~ 10~^ eV the situation is more complex since the resonance happens 
when Compton scattering is neither efficient nor inefficient. The evolution of the 
distortions has to be studied by numerically evolving the spectrum with the Kompaneets 
equation. We believe this is beyond the scope of this paper. 

Whenever Compton scattering is efficient, the photon distribution acquires, in 
short time scale after the resonance, a Bose-Einstein shape characterized by a chemical 




(35) 




(36) 



Fig. 3. 
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potential. For small distortions, this is given by [29] 

^'r = ^ " ^ ^^l'o(^0,nG) , (37) 

where Sp/p and 6n/n are the fractions of the CMB energy and photon number converted 
into ALPs during the resonance. Note that, in contrast to the case often studied in the 
literature, resonant photon-ALP conversion absorbs relatively more energy than photon 
number and therefore produces negative chemical potentials. 

Inverse Double Compton scattering (DC) and Bremsstrahlung (BS) occur much less 
frequent than Compton scattering but in contrast to the latter can change the photon 
number. Therefore, on a longer time scale they can absorb the necessary photons to 
recreate a pure Planck distribution, i.e. they can erase the chemical potential. For small 
distortions, the evolution of p can be approximated by [29] 

toe = 1.06 X 10^ ze^/^ s. ; tes = 3.73 x lO^^e^^^^ s , (39) 

where we have taken the values Yp = 0.25, Q^h"^ = 0.0223 for the ACDM model and 
defined zq = z/10^. During radiation domination we can write the expansion time as 
t = {2HY^ ~ 2.38 X lO^^e"^ s. We have computed for each ALP mass the time at which 
the resonance takes place and the evolution of the chemical potential. The FIRAS data 
sets the bound on the value of the chemical potential today yUo < 9 x 10~^, which gives 
the exclusion range denoted p in Fig. 3. 

The bounds from the distortions of the CMB spectrum vanish very fast for masses 
above ~ 0.1 meV because Double Compton scattering and Bremsstrahlung become very 
efficient. However, there is a further bound we can consider for resonances happening 
in the post BBN epoch. During ALP production a fraction of the energy stored in the 
CMB is transferred to ALPs, which immediately decouple from the thermal bath. If 
the ALPs produced in this way have sufficiently small masses to be relativistic during 
the epoch of the CMB formation (i.e. roughly m,^ <^ 1 eV) then they behave as a non- 
standard contribution to the radiation energy density. Therefore, the radiation energy 
density measured from CMB anisotropies would keep track of the ALPs contribution. 
A comparison of the radiation energy density during Big Bang Nucleosynthesis (not 
affected by ALPs resonant production) and CMB decoupling gives an upper limit on 
this contribution, namely x < 5p/ p = 0.2 [9, 30]. In our case this translates into a bound 

gioiBl^cV^' < 1.4 . (40) 

We show our combined constraints in Fig. 3. In this figure we have used the full 
expression in Eq. (35) since above ~ 0.1 meV the photon-ALP conversion probability 
bound cannot be considered small (also, for any mass a fraction of 0.2 of the photons 
converted into ALPs requires P^^tp of order 1). 
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Figure 3. Bounds on the ALP parameter space from distortions of the CMB 
blackbody spectrum caused by resonant j (j) oscillations. The ALP mass determines 
the time of the resonant transition. In the region labeled as "post-recombination" 
the transition happens in the post-recombination epoch whereas outside this region 
the transition occurs before recombination and CMB last scattering. However, in 
the region labeled as "weak-coupling" the transitions happen when the plasma is 
so weakly coupled that it cannot process the distortions. In the region labeled n, 
Compton scattering would restore a Bose-Einstein spectrum with a chemical potential 
that can be erased through inverse double Compton scattering and Bremsstrahhmg. 
For m0>O.l meV all distortions are erased but the produced ALPs still contribute 
to the cosmic radiation density and can be constrained by comparing the number of 
effective number of relativistic species at BBN and CMB decoupling. Here g is the 
ALP-photon coupling constant and {B^y^'^ an sky average of the comoving magnetic 
field during the resonance. 



7. Discussion and conclusions 

In this paper, we have calculated bounds on photon- ALP oscillations in the primordial 
magnetic field, deriving updated constraints from the high precision CMB spectrum 
data collected by the FIRAS instrument on board of COBE. A previous study [6] was 
derived in the pre-COBE era and it lacked a detailed treatment of the effects of the 
plasma medium on the photon- ALPs oscillations. This has motivated us to re-evaluate 
the bounds. We obtain limits on the product of the ALP-photon couphng g times the 
sky and polarization averaged magnetic field (B^)^/^, 

(?(fi2)i/2< iQ-n _ iQ-n GeV-^nG , (41) 

for ALP masses between 10^^"^ eV and 10""^ eV. Slightly weaker bounds were also derived 
for higher ALP masses. 

Our bound nicely connects with the one obtained in [11] for ALP masses less than 
10"^^ eV, considering only non-resonant conversions today. Our argument allows to 
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extend the sensitivity of CMB measurements also in a region of the parameter space that 
was previously unconstrained. We mention that ALPs can also be thermally produced 
in Early Universe via processes like c-f — > e0, forming a relic background. The (/>77 
vertex would also allow for radiative decays of these particles. In this case, additional 
cosmological bounds can be obtained for ^ eV, as studied in [31]. 

There are several interesting proposals to measure the CMB spectrum with higher 
sensitivity than the FIRAS instrument. The DIMES (Diffuse Microwave Emission 
Survey) proposal [32] aims to probe the region of smaller frequencies (2-100 GHz 
corresponding to 0.035 < ojq/Tq < 1.8) with sensitivity comparable to FIRAS. 
Unfortunately, as discussed above, the distortions grow with increasing frequency, 
so DIMES would not significantly strengthen bounds on axion-photon oscillations, 
in contrast to constraints on hidden photon - photon mixing for which conversion 
probabilities are inversely proportional to the photon energy. For this reason, the 
proposed FIRAS II [33] would be more relevant for photon-axion mixing constraints 
since it would shrink the FIRAS error bars by almost two order of magnitude and reach 
much higher frequencies (60-3600 GHz, or 1 < ujq/Tq < 65). If such a project is realised 
the sensitivity to photon-ALP oscillations will very likely improve by more than an order 
of magnitude in g{B'^Y^'^. A slight improvement on our bound in the pre-recombination 
era could be achieved thanks to new experiments, like ARCADE (Absolute Radiometer 
for Cosmology, Astrophysics, and Diffuse Emission) which would constrain the chemical 
potential of the CMB spectrum down to 2 x 10~^ [34]. 

Let us recall that recent results from the CAST experiment [35-37] give a 
direct experimental bound on the ALP-photon coupling of (7<8.8 x 10^^^ GeV^^ for 
m<^^0.02 eV, slightly stronger than the long-standing globular-cluster limit [38]. For 
ultra-light ALPs {m^ <, 10~^° eV) a stringent limit from the absence of 7-rays from 
SN 1987A gives ^ < 1 x 10"" GeV-^ [39] or even ^ < 3 x lO^^^ GeV-^ [40]. We stress that 
without direct evidence for a primordial magnetic field, our bounds on g{B'^)^^'^ do not 
allow to constrain directly the coupling constant g. However, if a primordial magnetic 
field would be found with values close to the current upper bound, the resulting CMB 
limit on g for ^ 10~^ eV would overcome the barrier placed by current experimental 
and astrophysical bounds. Conversely, if ALPs will be eventually discovered improving 
the current sensitivity of the solar axion helioscope CAST, or with new techniques [41- 
43] in laser experiments like photon regeneration [44-47] or laser polarization [43,48], 
our cosmological argument will provide a complementary constrain on the strength of 
the primordial magnetic field. 

Moreover, we note that our limit applies also to the so called chameleon-like scalar 
particles mixing with photons [49,50]. Such particles are introduced to explain the 
acceleration of the Universe, as inflation or dark energy fields [51], or to cause variations 
in the fundamental constants [52,53]. In general, chameleons have properties that 
depend on the environment. For this reason, they can evade astrophysical and CAST 
bounds, since in the dense stellar environment they become so heavy that they can 
not be produced in the usual reactions [18]. In these cases one has thus to rely on 
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the limits from the laboratory experiments, which take place essentially in vacuum. 
For a wide class of chameleon models, the laser PVLAS experiment would rule out 
values of x 10~^ GeV~^ [54] (see also [55-58]). Strong bounds also come from 
observations of starlight polarization [59]. In the expanding Universe, given the low 
plasma density, these particles would behave essentially as standard ALPs, so that our 
bounds can be directly applied to this case. In this sense, if primordial magnetic fields 
will be discovered, our argument would rule-out the recently proposed mechanism of 
chameleon-photon conversions to explain the observed supernovae la brightening [60]. 

Furthermore, let us point that the kind of analysis performed in this work can be 
also extended to other particles having two-photon vertices and sub-eV masses. For 
instance, a very similar analysis could be performed to massive spin-2 particles [61] like 
Kaluza-Klein gravitons [62]. 

As a final remark, our study is relevant for recent studies of conversions of high- 
energy gamma rays into ALPs in the intergalactic magnetic field in relation to the 
observed cosmic transparency of high energy gamma sources [63-65] or to ultra high- 
energy cosmic ray propagation [66]. If eventually a primordial magnetic field will be 
discovered, a definitive verdict on the impact of these fascinating mechanisms would 
require a combined analysis including our new CMB constraints. In this sense, it is 
nice to realize that our cosmological limits obtained with microwave photons could 
have relevant consequences for signatures of ALPs in high-energy gamma sources. This 
confirms once more the broad range of energies potentially accessible through axion-like 
particle searches. 
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